function y = edengsch(x)
% EDENGSCH function
% x0=[0,0,...0];
n=length(x); y=16;
for i=1:n-1
    y=y+(x(i)-2)^4+(x(i)*x(i+1)-2*x(i+1))^2+(x(i+1)+1)^2;
end


